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We study a possible gravitational vacuum-effect, in which vacuum-energy variation is due 
to variation of gravitational field, vacuum state gains gravitational energy and releases it 
by spontaneous photon emissions. Based on the path-integral representation, we present a 
, general formulation of vacuum transition matrix and energy-momentum tensor of a quantum 

■ scalar field theory in curved spacetime. Using analytical continuation of dimensionality of 



the phase space, we calculate the difference of vacuum-energy densities in the presence and 
absence of gravitational field. Using the dynamical equation of gravitational collapse, we 
compute the rate of vacuum state gaining gravitational energy. Computing the transition 
amplitude from initial vacuum state to final vacuum state in gravitational collapsing process, 



^SJ ' we show the rate and spectrum of spontaneous photon emissions for releasing gravitational 

energy. The possible connection of our study to the genuine origin of gamma ray bursts is 



discussed. We compare our idea with the Schwinger idea for Sonoluminiescence and contrast 



. our scenario with the Hawking effect. 

O ; PACS numbers: 04.62. -fv, 04.70.Dy 
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^ , I. INTRODUCTION 

> : 
•1—1 , 

^ . The issue of quantum field theories of elementary particles in curved spacetime has played a 

Cd \ tremendously important role in understanding quantum phenomenon of particle creations from the 

vacuum, when gravitational field is present. In quantum field theories, the vacuum state consists of 
a large number of virtual particles, that are quantum-field fluctuations whose energy-momentum 
are off-shell. The vacuum energy (zero-point energy) is attributed to the energy-momentum of 
virtual particles. These virtual particles interact with a external gravitational field via their energy- 
momentum. Such interaction results in vacuum-energy variations and particle creations from the 
vacuum. 

The Hawking radiation is the phenomenon of particle creations occurring around black 
hole's horizon. The quantum radiation of the Hawking type is rather general in curved spacetimeP]. 
Although, such quantum radiation is too small to be detected in the present Universe, its impact 



'Electronic address: IxueQicra.itI 



2 



on theoretical understanding of particle creations due to a static gravitational field is indeed far 
reaching. On the other hand, it is important to study the phenomenon of particle creations in a 
non-static curved spacetime. In fact, the rate and spectrum of particle and antiparticle creations in 
extremely early Universe are possibly related to the CMB pattern and the large-scale structure of 
the present Universe, which are most exciting arena of theoretical and observational physics today. 
Beside, in a gravitational collapse process, violent variation of gravitational field possibly causes a 
large number of particle creations, which might be account for most energetical events of gamma 
ray bursts 4]. This is the topic that we attempt to further discuss in this article. 

The phenomenon of particle creations from the vacuum was first studied by Euler, Heisenberg 
and Schwinger in the quantum field theory of electromagnetic dynamics (QED) |5|. A strong 
external electric field, its strength is larger than m^c^/h, greatly reduces the energy-mass gap 
(~ 2nie) of charged virtual particles in the QED vacuum, so that virtual particles undergo a 
quantum tunneling process, leading to large pair-productions of electrons and positrons. The 
mechanism and phenomenon of electron and positron pair-productions have advocated a numerous 
studies both in experimental physics ^ and theoretical physics^?], as well as in astrophysics 8j • In 
addition, it is worthwhile to mention that an external magnetic field can induce vacuum decay 
leading to spontaneous photon emissionsy]. The dynamics of this effect, which is completely 
different from the Schwinger one, is that the external magnetic field modifies the vacuum-energy 
spectrum so that the vacuum state gains magnetic energy and becomes unstable. This effect can 
be possibly tested in laboratories and account for the phenomenon of anomalous x-ray emission 
from pulsar s IkI . 

Casimir IjJ first considered that boundary conditions (two conducting plates) modify the energy 
spectrum of the vacuum state and vacuum-energy variation is the Casimir energy 5<?|casimir- An 
attractive force between two plates is observed as the Casimir effect The reasons for the 
Casimir effect are that the vacuum state gains the Casimir energy and becomes energetically 
unstable, quantum-field fluctuations result in releasing the Casimir energy. Recently, there is 
much theoretical and experimental attention on the dynamical Casimir effect [isl. In such an 
effect, boundary conditions are dynamically time-dependent, the Casimir energy is released by 
spontaneous photon emissions. 

Sonoluminiescence [l^. Iisj l is another most interesting vacuum effect which shows a flash of 
spontaneous photon emissions, when gas bubbles in water collapse, driven by a sound-wave. The 
vacuum-energy variation is due to the variation of dielectric constant, rather than the modification 
of boundary conditions as in the Casimir effect. We will have a more detailed discussion in this 



article Q. 

The mystery of energetic sources generating gamma ray bursts [ifll is a prompt emission 
(seconds) of extremely huge energy output (~ lO^^ergs for isotropic emission) from rather compact 
sources (~ lO^cm) at cosmological distance (z ^ 0(1))- These have stimulated many studies in 
connection with electromagnetic properties of black holesj^, 3)^^- Various astrophysical scenarios 
are discussed in literatures Il9|l. It is worthwhile to mention that via the Schwinger mechanism of 
electron-positron pair productions, the "dyadosphere" [3, 120|] of photons, electrons and positrons, is 
formed during the process of gravitational collapse of a massive star with electromagnetic structure. 

In this article, we want to study a possible gravitational vacuum-effect, in which vacuum-energy 
variation is due to the variation of gravitational field, vacuum state gains gravitational energy and 
releases it by spontaneous photon emissions. In section (JJ), using path-integral representation, 
we present a general formulation of quantum scalar field theories in curved spacetime, where the 
vacuum states are defined, the transition matrix from vacuum to vacuum, vacuum-energy spectrum 
and vacuum energy-momentum tensor are obtained. In sections pil|) and HIV|) . we specify two 
static observers respectively in the Schwarzschild and flat spacetime; we analyze the eigenvalues 
(spectrum) of the transition matrix and define the phase space of vacuum states; by using analytical 
continuation of dimensionality of the phase space, we obtain vacuum-energy density. In sections 
(fVj) and ()VI|) . we calculate vacuum-energy density and discuss why the characteristic energy-scale 
should be the ultraviolet cutoff for considering the difference of vacuum-energy densities in the 
presence and absence of gravitational field. In section ()VII|) . we define vacuum energy with respect 
to two static observers and compute the difference of vacuum energies in the presence and absence 
of gravitational field, which shows vacuum state gains gravitational energy. In section ()VIII|) . we 
adopt a simple model describing the gravitational collapse of a massive shell that is infinitesimally 
thin, and using the dynamical equation of gravitational collapse, we compute the rate of vacuum 
state gaining gravitational energy. In section ()IX|) . based on the transition amplitude from initial 
vacuum state to final vacuum state at each step of gravitational collapsing process, we compute the 
rate and spectrum of spontaneous photon emissions for releasing gravitational energy. In sections 
(pCj) and ()XI|) . we compare our proposal for gamma ray bursts with the Schwinger proposal for 
Sonoluminiescence and contrast our scenario with the Hawking effect. In the final section (|XII|) . 
we make some remarks on this preliminary study and the possible connection of our proposal to 
the genuine origin of gamma ray bursts. 
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II. GENERAL FORMULATION. 

We assume that the structure of spacetime is described by the pseudo-Riemannian metric g^y 
associated with the line element 

ds^ = g^vdx^dx^ , ;U, = 0, 1, 2, 3, (1) 

and the spacetime point is coordinated by x = (x*^, x*) = {t, x). The special geometrical symmetries 
of the spacetime S are described by using Killing vectors i^'^, which are solutions of Killing's equation 

C!.g^u{x) = Q, i^;u + iu;t, = ^, (2) 

where is the Lie derivative along the vector field orthogonal to the spacelike hypersurface 
S( {t =constant) of the spacetime S. A static observer O is at rest in this hypersurface S^. 

We consider that at the initial time (tin = —^t/'^), the spacetime is asymptotically flat, described 
by asymptotically flat geometry g^i^ in Eq.(^; while at the final time (tout = +St/2), the spacetime 
is curved and stationary, described by a non-trivial geometry g^i, in Eq.(^. The characteristic 
time-scale 5t of such variation of spacetime geometry is supposed to be much larger than the 
characteristic time-scale (e.g., l/nie) of a quantum-field transition, i.e., 5t ^ 1 in the unit of 
the quantum time-scale. We attempt to study quantum-field fluctuations interacting with the 
variation of spacetime geometry. A specific model for such variation of spacetime geometry and 
quantum-field transition will be presented in due course. 

In order to clearly illustrate the physics content, we first consider a complex scalar field (f) in 
curved spacetime. The simplest coordinate-invariant action is given hy {h = c = G = k = 1) 

S = l f d^xV^U'cP^^cP* + (m^ + , (3) 



2 

where m is particle mass and TZ the Riemann scalar. The quantum scalar field (j) can be in principle 
expressed in terms of a complete and orthogonal basis of quantum-field states Uk{x): 



(/)(x) = (akUk{x) + alul{x)j , a^, 4, = Sk,k' (4) 



where a\, and are creation and annihilation operators of the k-th quantum- field state Uk{x). 
This quantum field state obeys the following equation of motion, 

(A^ + m2 + e7^)nfc(x) =0, (5) 
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and appropriate boundary conditions for selected values of k. In Eq.©, A^. is the Laplacian 
operator in curved spacetime: 

1 



-X - -^d^iV^g'^'d,]. (6) 

In Eq.@, we assume that Ufc(x) are positive energy (w) states, satisfying 

C^Uk{x) = —ioJUkix), u; > 0, (7) 

with respect to the timelike Killing vector field ^'^ © associated to the static observer O rest in 
hypersurface St [t =constant) of the spacetime S. 

At the initial time (tin = —Stf^), the spacetime is approximately flat, quantum-field states 
Uk{x) © are asymptotically free states Uk{x), obeying Eq.© for the asymptotically free geometry 
Qf^i^. Then, the quantum scalar field (f){x) is an asymptotically free field in the hypersurface '^-st/2 
of the asymptotically free spacetime S: 

4>m{x) = (akUkix) + aluKxfj , a^, 4' = ^k,k' (8) 

k 

where x G S_5t/2 5 are creation and annihilation operators of the A;-th asymptotically free 

quantum- field state Uk{x). Corresponding Lie derivative along the Killing vector (j^J is dt, positive 
energy states are Uk{x), satisfying Eq.Q. Then we may construct the standard Minkowski space 
quantum vacuum state |0,in): 

afc|0,in) = 0, (0,in|4 = 0. (9) 

|0, in) is an initial quantum vacuum state at the initial time (tin = —St/2) with respect to the static 
observer O rest in hypersurface S_5j/2 of the spacetime S. 

At the final time (tout = +St/2), the spacetime S is curved, described by a stationary non- 
trivial geometry Qf^^. We assume that quantum-field states Uk{x) © are asymptotical states 
Ukix), obeying Eq.© for the stationary geometry g^^. In the hypersurface ^+st/2 of the spacetime 
S, the asymptotical quantum scalar field 0out is expressed in terms of Uk{x): 



0out(2;) = ^ [dkUk{x) + a\ul{x)j , dk, a\, = 6k,k' (10) 



k 



where x G '^+5t/2j o,k ^^'^ are creation and annihilation operators of the A;-th quantum- field 
state Uk{x). Corresponding Lie derivative along the Killing vector is ((21), positive energy states 
are Uk{x) satisfying Eq.Q. Then we may construct the quantum vacuum state |C),out): 



afc|0,out) = 0, (0,out|4 = 0, (11) 
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in curved spacetime. |0,out) is an final quantum vacuum state at the final time (tout = +5t/2) 
with respect to the same static observer O rest in hypersurface '^+st/2 of the spacetime S. 

It is worthwhile to note that </>out(2;)({'Ufc(x)}) are not asymptotically free field(states), instead 
they are asymptotical field(states) in the presence of external and stationary gravitational field, so 
that the final vacuum state |0,out) is different from the initial vacuum state |0, in) Q. Such 
a difference is not only a unitary phase. The final vacuum state |0,out) may not necessarily 
be measured as devoid of particles, in contrast to the initial vacuum state |0,in) defined by Eq.Q 
relating to the asymptotically free field (j)\-a{x). In fact, as will be shown, the final vacuum state 
|0, out) is a quantum-field state of particle and antiparticle creations upon the initial vacuum state 
|0,in). This indicates gravitational field interacting with quantum-field fluctuations of positive 
and negative energy states of the initial vacuum |0,in), and the quantum scalar fleld evolves 
throughout intermediate quantum-field states (j){x,t) Q) for —5t/2 < t < +5t/2. We speculate 
that this evolution is adiabatic for 5t being much larger than quantum-field transition time. 

To deal with all possible intermediate states, represented by (j){x) or Uk{x) in Eq.@ for —6t/2 < 
t < +6t/2, we use path-integral representation to study the transition amplitude between the initial 
vacuum state and final vacuum state: 

(O,out|0,in) = J[V(l)V^*]exp{iS), (12) 

where 

J [VW] = n_5i/2<i<+5i/2n,-es, / tW{x, t)]. (13) 

The intermediate quantum-field states contributions to the transition amplitude (fT^ can be for- 
mally path-integrated, 

(O,out|0,in) = det"^ (7W) , 7W = + + ^7^. (14) 

This result clearly depends on the initial vacuum state (j)i-a (jHl) and final vacuum state 0out 
which are not explicitly written. The effective action S^s is defined as 

^eff = — i ln(0, out|0, in), (15) 

which relates to the phase of the S-matrix transition from the initial vacuum state |0, in) to the 
final vacuum state |0,out). The averaged energy-momentum tensor {T^jj) of the quantum-field 
vacuum is given by: 

(r^^(x)) = ^T^. (16) 
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Main efforts are calculations of the effective action ()15() and energy-momentum tensor (|16() in this 
article. 

In order to evaluate the path- integral (|12|) over all intermediate quantum- field states (^(x) 
it is convenient to introduce operators and defined on the states \x) and \k): 

X^\x)=x^\x)- k^\k) = k^\k). (17) 

They enjoy the canonical communication, 

[X^,k,] = -ig^,. (18) 

The states \x) and \k) satisfy: 

{x\x) = 5{x — x'), j dx\x){x\=l 

{k\k') = 27r6{k-k'), J dk\k){k\ = 1, (19) 

and intermediate quantum- field state Uk{x) can be represented as 

Uk{x) = {x\k). (20) 

Using these matrix notations, we write the operator M{X,K) (|14j) as a hermitian matrix 

= {k\x){x\M{X,k)\x'){x'\k') 
= u%{x)M.{x^ K)uk'{x')6{x — x'), 

= ul{x)M{x, K)uki{x')5kk'5{x - x), (21) 

where 

{x\M{X, k)\x') = M{x, k)6{x - x'), (22) 

is diagonal in the coordinate space {x}. In the representation of {uk{x)}, the hermitian matrix 
■Mk'ki is also diagonal in the energy- momentum space {k}, since {uk{x)} are eigenstates of the 
operator M ©. By normalizing the quantum field (p, we define the normalized diagonal element 
of the matrix (|21jl as 

Xl{x) ^ ^^ul{x)M{x,k)uk{x), (23) 
and formally compute the effective action Ses given in Eas. (|15ll4p : 



iSeS = — ln(0, out|0, in)z5eff 
-d'^xd'^k 

-d'^xd'^k 



/. — d xd k 2/ M 



5^^^^1nAi(x)|inj , (24) 
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where the {A^jout and {A^}in are the diagonal elements (^5]) . The {A|}out is in terms of the final 
vacuum state and geometry g^^y^ whereas the {A|}in is in terms of the initial vacuum state c/^in and 
geometry g^y. 

Allout = |~ ] .,o Ul{x)M{x,k)uk{x), (25) 
Ailin = |_ , ul{x)M{x,k)uk{x). (26) 

The operator A|,(x) (PS)) and the number of quantum-field states / y/—gd^xd^k/{2'KY are invariant 
in arbitrary coordinate systems, later is the Liouville theorem for the phase-space invariance. 

By using Eqs. (|l(i|) and (|24|) . the variation (T^i/(x))difr of averaged energy-momentum tensor in 
the time interval 5t is given by the difference: 

{T^''{x))Am = {T^''{x))out - {T^''{x))in, (27) 

where (T'^'^(a;))out is the averaged energy-momentum tensor computed by the variation of geometry 
9fiu ifT^ - corresponding to the final quantum vacuum state (pout, whereas {T^^{x))\yi is the averaged 
energy-momentum tensors computed by the variation of geometry g^j^y (|T6|) . corresponding to the 
initial quantum vacuum state (pla- 
in general, using the definition (|24|) . we formally calculate the averaged energy-momentum 
tensor (Tf^^^x)) as, 

(r^^(x)) = (r^'^(x))(^) + (r^''(x))(2), (28) 



where 



(r^^(x))(^) = ig^'^ix) J ^ln(Al(x)) (29) 



(27r) 

(T-(x))(^) = 2^ I ^^^^ 



(30) 



In these equations 



'-ipB = \^)r{x)6\x-y)- 



d^x^5\x-y)f{x) = f{y), (31) 



for an arbitrary function f{x). The calculations of energy- momentum tensor ()29I3U|) are main tasks 
in the following two sections. 
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III. VACUUM ENERGY-SPECTRUM. 

In order to illustrate physical idea in a mathematically tractable way, we model a collapsing mas- 
sive star as an infinitesimally thin and spherical shell. This massive shell separates the spacetime 
into two regions: (i) internal region S described by the flat geometry g^^, = (1, —1, —1, —1), 

ds^ = dtl-dr^-r^dn, r < R; (32) 

and (ii) the external spacetime S described by the stationary Schwarzschild geometry, 

ds'^ = gttdt^ +grrdr^ — r^dVL^ r > R 

9tt = -(grr)'^ = g = il-'^), (33) 

where r, 9, (j) are spherical-polar coordinates, dO, = dO'^ + sin'^ Odcj)^ , t and are the Schwarzschild- 
like coordinates in the external and internal spacetime respectively. R indicates the radial position 
of the shell. M is the total mass-energy of the shell. 

At the initial time tin = —St/2, the shell radius is R and a static observer O is located at 
x{R — 0^, fi) G S_5(/2 ill the internal spacetime S. His four velocity and Killing vector are 
given by, 

= (1,0,0,0), ^^ = (1,0,0,0). (34) 

The initial quantum field is (pi^ Q and quantum vacuum state is |0,in) Q. After St, the shell 
gravitationally collapses and its radial position R moves inwards to i? — SR. At the final time 
*out = +St/2, the same static observer O turns out to be in the external spacetime S. His four 
velocity and Killing vector are then given by, 

u^ = {9T{r), 0,0,0), C^. = {9tt{r), 0,0,0). (35) 

The final quantum field is (j)o^t (|10|) and quantum vacuum state is |0,out) (jllj) . In the following, 
we respectively compute the vacuum-energy spectrum of the initial and final vacuum states. 

In the external spacetime, the Riemann scalar TZ = and the Laplacian operator Q is given 

by: 

^ ^ dt^ dr ^ dr 

^ dt'^ ^ I dr"^ r dr j 
- + 7^' (36) 
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where L?' is the angular momentum operator, g** = [gtt] ^ and g"^^ = (grr) ^■ 
The appropriate basis of quantum field states (jlOj) is chosen as 

Uk{x) = {t,r,e,^\u;,kr,l,m) = RUr)Yi„X0,4>)e"''\ (37) 

where Yim{6, (p) is the standard spherical harmonic function: L'^Yim{0, (p) = l{l + l)Yim{G, (p) and k 
indicates a set of quantum numbers (w, kr, I, m). uj is the energy-spectrum and the radial momen- 
tum kr will be defined soon. From Eq.©, the radial function Riui{r) obeys the following differential 
equation for r > R > 2M, 

2M, 



g^W + g'^^P,+i=^kr-Vi 



RUr) = 0, (38) 



where the hermitian radial momentum operator 

ir dr ^ ^ \ dr"^ r dr 



kr = -^r, fc2 ^ - + - , 39 



and potential 



^^ = ^7^ + — + "^'- (40) 



Ea. ()38() is exactly equivalent to the Regge and Wheeler equation. The radial function Rioj{r) 
is an orthogonal and complete basis, asymptotically behaves as the Hankel function hii^{r) for 
r ^ R > 2M. The matrix operator M in Eq. (|2Hj) is given by, 

M{r, K) = - 9'''kr - ^^kr + Vi. (41) 

Eas. ()37l41|) define a complex eigen- value problem to find the energy spectrum of the final vacuum 
state in an external gravitation field. The imaginary part of the operator M (r, kr) (|41j) results in 
the quantum radiation of Hawking type in curved spacetime, as discussed in ref.j^. 
The diagonal matrix A| (^H)) is given by: 

Ai(x)|out = -g'^oo' - f'kl + Vi- i'^kr, (42) 
where we define the values of "radial momentum" kr and k^. of the quantum field state Ukix): 

_ ulkrUk _ R^jiKRuA^ ,.„s 

2 _ ulkrUk _ Rljik^Ruii , . 

In the internal spacetime 5, the Laplacian operator ® is given by: 
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The appropriate basis of quantum field states © is chosen as 

Uk{x) = {to,r,9,cl)\LOo,kro,l,m) = Ri^,{r)Yim{0,<P)e-'^''' ■ (46) 
The radial function RiujQ{r) obeys the following differential equation, 

Ri^oir) = 0. (47) 



^ m 



The radial function Riujq{t) forms a orthogonal and complete basis. The matrix operator M in 
Ea. (PH|) is given by, 

M{t, kr) = ^ + k^r + ^^4^ + m^- (48) 
otq r 

The diagonal matrix (|23j) is given by: 

Ai(x)|in = -ojI + k% + ^^^^ + m^, (49) 
where we define the values of "radial momentum" kj-o and /c^q of the quantum field state Uk{x): 

^ 1^ ~ \R.,l\' ' ^^^^ 

U2 _ Ul^r^k _ RlolkyR mol 



'^rO — |-_ |2 — IP 12 • 



For m = in Eq.(refeq2f), R^^^^i = 2a;oJz(wor), ji{u)or) is the spherical Bessel function. As a 
particular case, we adopt the spherically symmetric solution / = and m = in the differential 
equation Ea. H47() . such that the spherically symmetric solution -Roa;o(^) ~ e*^'''"/r is the eigenstate 
of the radial momentum operators k^ and k^. (IHU]) . The eigenvalues kro (jSn]) and /c^q related 
to radial momentum of the spherical quantum field (/ = 0) . This indicates that the values of "radial 
momentum" k-r (|43() and k"^ (|44|) are consistent with the radial momentum of quantum field states. 

IV. VACUUM ENERGY DENSITY 

Armed with the energy- momentum tensor {T^'^{x)) ()28() and the energy-spectrum A|(2;)|out (02) 
in the external spacetime S, we are ready to compute the energy-density, 

(r/)out = 5«(T")out. (52) 

Given the volume of spherical shell Airr'^drdt, we have the number of quantum field states in the 

energy-momentum phase-space, 

d'^k f d'^k± f dwdkr 1 [ dujdkr 

2^ / 7^' (^2) 



(27r)4 J {2^fJ (27r)2 47rr2 ^ 7 (27r) 
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where kr is the "radial momentum" defined in Ea. ()43() and k± are the transverse momenta, per- 
pendicular to the radial direction. 

Starting with the first part {T^'^ (x))^^^ of the energy-momentum tensor p9|). we compute the 
vacuum energy-density (fS^ as: 



47rr2 



E 

Lm 



dwdkr 



ln(A2). 



Using the identity: 



we are able to write the vacuum energy-density ((54|) . 



In j"^ Lis{b+ie) _ ^is{a+ie)\ _ 

b Jo s y J ' 



(54) 



(55) 



E 

l,m 
1). 



diodkj. 



°° ds 
s 



(56) 



The second term indicated by (A| ^ 1) is the same as the first term with A| ^ 1, and this term 
is a constant. The logarithmatic function in Ea. H54|) is represented by an s-integration in Ea. H56|l 
and infrared convergence at s — > is insured by ie prescription (e — > 0). 

Introducing the variable ^ = —ito (the Wick rotation) and the integral representation, 



(2^)' 



1 



(57) 



for Im(/3) < 0, we express Ea. (|^ as, 



Sttv 



E 

Lm 



+ (Ai - 1) 



(27r) ^^i^ Jo s3/2 



°° ds is(Xl+ie) 



(58) 



where and henceforth, is given by Ea. H42() without the cj^-term. 

In order to compute the integration over "s" in Ea. (|58() . we introduce a complex variable 
z = —1/2 -|- 6 {\6\ 0), and use the following integral representation of the r(z)-function by an 
analytical continuation for Im(a) > 0: 







ds = {-iayT{z). 



(59) 



This analytical continuation is equivalent to analytical continuation of dimensionality of the 
momentum-space / in Ea. ()57() . In the neighborhood of singularity, where \6\ — > and 
z —1/2 in Eq.((Sni), we have 



T{z) = —2^/tt, a 



I gk"^ + Vi- i^^kr. 



(60) 
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Substituting Eas. H59l60|) into Ea. (|5S|) . we obtain the vacuum energy-density: 

1/2 _ r dkr ' 



47rr2 



E 



{2n) 



I gk^ + Vi-i^^kr, 



(61) 



where the "radial momentum" integration / dkr integrates from "0" to a ultra violate cutoff A. 

Now we turn to the computations of the second part (T^'^(x))(^) of the energy-momentum tensor 
pUj) . Using Eas. H30l52|) and exchanging momentum- integration and metric- variation, we write 



't\{2) 



2gtt 



^gtt{x) 



d^k 

(2^)4 



Hxlix)) 



(62) 



where the bracket [• • •] was computed, as shown in Ea.((^. In addition, we have S{g*'^gu) = 0, 
g^^Sgtt = -gttSg^^ and 



9tt- 



^9tt \/g^ <55** Vi" 

Using Ea. (j55)) and relationship we take metric-derivative in Ea. (|H^ and obtain the result: 



(63) 



(T,*)(^) = (T*)«. 



Thus the vacuum-energy density in the external spacetime S is 



(64) 
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1/2 



E 



I gk'^^+Vi-i'^kr 



l,m ^ ' ' 

This is the vacuum energy-density corresponding to the final vacuum state Bout- 
in the internal flat spacetime 5, the vacuum-energy density H65() is reduced to. 



(65) 



47]-r2 



(2vr) 



(66) 



This is the vacuum energy-density corresponding to the final vacuum state (\)-xa. We adopt the 
transverse momenta A:_l to replace the quantum numbers (/, m) of angular momenta: 

/(/ + 1) /• d^k^ 1 



47]-r2 



E- 

Lm 



(67) 



Eg . (1661) can be written as 



d^k 
(2^ 



(68) 



where k"^ = k^Q + d^k = dkrod'^k± and integration range [— A, A]^. Eg . (1681) is consistent with 

the vacuum-energy density for two components of quantum scalar field. 
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Comparing the vacuum-energy density (T/)out (|fc)5(l in the Schwarzschild geometry with the 
vacuum-energy density {Tl)i^ (|66jl in the flat geometry, we find the vacuum-energy of the quantum- 
field fiuctuations (virtual particles) of the vacuum couples to the gravitational field in the following 
aspects: (i) an addition term 2M/r^ in the potential Vi 1)40(1 and an imaginary term i2Mkr /r'^\ (ii) 
the "radial momentum" kr and k"^ are modified from Eqs. ()4l'{l44|) to Eqs. ()5()l51j) : (iii) modification 
gkl k^Q. This shows the radial modes described by kr directly interact with gravitational field. 
While, the transverse momenta k± in Ea.H65|) are the same as its counterpart in Ea.H66|). showing 
the transverse modes described by k± do not directly interact with gravitational field. In following 
sections, we attempt to calculate the difference of vacuum-energy densities, 

(r*(x))difr = (r/(x))out - (69) 

where (T/(x))out is in terms of the final quantum vacuum state 0out and geometry g^i,; whereas 
(T/(x))in is in terms of the final quantum vacuum state (pi^ and geometry 5^;^. 

V. THREE ENERGY SCALES A, m AND r ^ 

The vacuum-energy densities Eas. (|65() and H66|) depend on the three very different energy scales: 
A, m, r^^ ( A ^ m ^ r^^). The scale M is considered the same as the scale r~^. One would 
expect that final results of vacuum-energy densities (T/)out and {T^)\a should be in terms 
of finite terms: r~'^,'mr~^'' ^m?r^'^ and m^, since all "divergent terms" containing A should be 
removed away. In what follows, we argue that A is a physical cutoff of its own right for calculating 
the difference (|69|) of vacuum energies in the presence and absence of gravitational field; "divergent 
terms" containing A should not be simply removed away. 

We know that a quantum field theory in the fiat spacetime, although its characteristic energy- 
scale is the measured mass m of particles, has involved high-energy modes in the computations of 
virtual particle contributions (loops in Feynman diagrams) to quantum corrections. The ultraviolet 
divergent terms arise from the contributions of these high-energy virtual particles. A ultraviolet 
cutoff A or another method is introduced to regularize these ultraviolet divergent terms, which are 
then consistently removed away by renormalizing quantum fields and parameters of theory. Only 
relevant contributions of virtual particles, whose wave-lengths are the order of m^^, are taken into 
account. This is due to the fact that high-energy virtual particles characterized by the cutoff A 
must not contribute to a physical process characterized by the energy scale m, that is much smaller 
than A. In a sensible quantum field theory, theoretical results corresponding to experimental low- 
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energy measurements, must be independent of the ultraviolet cutoff. The ultraviolet divergent 
terms can be consistently removed, if and only if the number of renormalized fields and parameters 
(for instance, coupling and mass) is exactly equal to the number of different types of ultraviolet 
divergent terms. This is normally guaranteed by internal symmetries of a renormalizable quantum 
field theory. 

The vacuum energy of a quantum field theory is a "divergent constant", since it is resulted 
from virtual particles whose energy range from to the cutoff A. In the flat spacetime, only the 
energy-difference of quantum states can be measured, as a consequence the vacuum energy can be 
discarded, by normal ordering of quantum fields. However, in studying the Casimir effect, we need 
to calculate the difference of vacuum energies in presence and absence of two conducting plates, 
since the presence of two conducting plates modifies vacuum state and its energy. This difference 
of vacuum energies turns out to be a measurable Casimir effect. The resultant Casimir energy is 
characterized by the size of the distance L between two conducting plates. This is because: 

• (i) those virtual particles of wavelengths comparable with L strongly impacted by two plates 
and their contributions to vacuum energy are modified, giving rise to the Casimir effect; 

• (ii) whereas those virtual particles of wavelength incomparable (either much smaller or 
greater than) with L are not much affected by two plates, so that their contributions to 
vacuum energy are same, and as a result these virtual particles have no contributions to the 
Casimir effect. 

The Hawking effect is certainly due to virtual particles in the vacuum interacting with an 
external static gravitational field near the horizon of a black hole. The reasons that this effect is 
characterized by the energy scale T ^ 1/M, which is the gravitational potential at the size of a 
black hole, are the following: 

• (1) the dynamics of such an effect is virtual particles quantum-mechanically tunneling 
through the gravitational potential near the horizon of a black hole; 

• (2) those virtual particles, whose energy is comparable with the gravitational potential 
around the horizon, have a large probability of undergoing the tunneling process, since 
there are more crossing energy-levels and virtual particles at these energy-levels can tunnel 
out of the gravitational potential into infinity; 

• (3) those virtual particles, whose energy is much smaller than the gravitational potential, 
have a small probability of undergoing the tunneling process, since there are not crossing 
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energy-levels and virtual particles at these energy-levels cannot tunnel out of the gravitational 
potential into infinity; 

• (4) those virtual particles, whose energy is much larger than the gravitational potential, are 
not disturbed by the gravitational potential and remain as virtual particles in the vacuum. 

Analogously, the quantum radiation in curved spacetime discussed in ref.j^] has the characteristic 
energy scale ~ 1/r, where is the gravitational potential that virtual particles tunnel through. 

We turn back to our case: we study a quantum field theory in an external gravitational field 
of a collapsing massive shell and compute the difference H69() of vacuum energies in the presence 
and absence of gravitational field. Gravitational field is a classical field of long wavelengths 0(r). 
In curved spacetime, the vacuum energy cannot be simply be discarded, since gravitational field 
couples to the vacuum energy. As we can see the vacuum to vacuum transition matrix Ea. H42() . 
the term 

shows that gravitational field couples to the energy-momentum {uj,kr) of virtual particles in the 
vacuum. High-energy modes of virtual particles have much stronger interactions with gravitational 
field than low-energy modes of virtual particles do. The cutoff A should be a real physical cutoff, 
determined by an energy scale, where the difference (|69j) of vacuum energies in the presence and 
absence of gravitational field vanishes. We speculate that the cutoff A seems to be the Planck energy 
scale Ap for the reasons: (i) virtual particles of high-energy up the Planck energy scale interact 
with gravitational field and (ii) we do not see any intermediate energy scale, where the difference 
H69|) of vacuum energies in the presence and absence of gravitational field vanishes. However, we 
leave the ultraviolet cutoff A as a parameter determined by the phenomenon of gamma ray bursts. 
We will also have a similar discussion on this point, in comparison with the ultraviolet cutoff K 
introduced by Schwinger for the phenomenon of Sonoluminiescence in section (pCj) . 

To end this section, we note two points: (i) the difference (|(i9|) of vacuum-energy densities 
(r/)out and (T/) in receives dominate contributions from high-energy modes of virtual particles, in 
particular, those modes at the ultra violate cutoff A; (ii) the back action of vacuum-energy variation 
to gravitational field is not considered. 
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VI. COMPUTATIONS OF VACUUM ENERGY DENSITIES 

We are interested in computing vacuum-energy densities Eq (|65|) in the presence of gravitation 
field and Eg . (1661) in the absence of gravitation field. Our aim is to find the difference (|69j) of 
these vacuum-energies. Therefore the modes of virtual particles that do not interact with the 
gravitational field should be discarded. As we can see from the potential term Vi H40|) in Ea. (|65() . 
the radial modes described by kr and the transverse modes by angular quantum numbers (l, m) or 
transverse momenta kx_ ()67() play very different roles in contributing to the vacuum-energy density. 
The radial modes of virtual particles of high energy up to the cutoff scale A, directly interacting 
with gravitational field, are most relevant, since they make dominate contributions to the difference 
of vacuum energies. While, the transverse modes do not directly interact with gravitational field, 
since the term k'j^ = l{l + l)/r^ in H65() is the same as that in ()66() . Thus, we separate the S-wave 
{I = 0) contributions from non S-wave {I ^ 0) contributions to energy-momentum tensors (|H5pH6j) in 
the following computations. In this section, we compute the S-wave {I = 0) contributions {TlY-~^ 
and {TiY^. The calculations of non S-wave {I / 0) contributions (T/)!^'' and {Ti)^^ will be 
presented in appendix A. 

We first compute vacuum-energy density (T/)in (|66|) in the absence of gravitational field. Sep- 
arating the S-wave contributions from non S-wave contributions we write Ea.H66|l 
as, 

(r/>in = (T/)(r + (T*)|f , (71) 

where 
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In Eg .(1731). we use the dimensional regularization of / dkro — > / S^'^^'>kro in Ea. (|72|) and analytic 
continuation representation: 

/ d^'^^^Ko^k^ + A-^ = r(-i)[^2](i+^), (74) 

where e — > and r(— e/2) ~ — 2/e. The divergent term 2/e represents the ultraviolet cutoff A of 
"radial momentum" ky-Q. 

We turn to compute the vacuum-energy density (T/)out in the presence of gravitational 
field. Separating the S-wave contributions from non S-wave contributions in Ea. H65|) . we have, 

(T/)out = (T*>i:° + (r*)it?, (75) 
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where 

{Tiy=^ = 2f^J ^^^9k?+^=o-^^-^K, (76) 
where Vi=o is Vi{r) (|40j) for / = 0. Changing the integrating variable in Eq. (j76j) : 

k-p ^ ^ Q ^ kf ^ ^ ky^ — kj, z — , (TT) 

we write Eas.(|7B|) as, 



(^/)iu? = ^/ ^^^k^ + m^ + Q^ (78) 



+ (79) 



where 



= — + — • 80 

In Ea.((7^. we use the formula ((7^ to compute the fc^-integrations in (T'/)^^^ ((75|l. analogously to 
the calculations of (T/)|=o (O- 

Using (T/)j=° dZSl) and (T/)[,=t° ((711) for e ^ up to 0(2/e), we find that the term is canceled 
in the difference {T^Y^^ dSH). To compute the kro- and fer-integrations in (r/)|=° (O and (r/)[,=t° 
(|78|) . we use the formula, 

dxVax^ + 62 = A^aA2 + 62 (81) 
JO 2 

&2 ^ /Aai/2 + VaA2+62\ 
+ b )■ 

For A ^ m ^ r^^, to the leading term (0(A2)) we approximately have, 

<^'>t° - 4^ (^) ^ («^' 

The vacuum-energy densities Eq.ljHH) and Ea. H83p are mainly contributed from virtual particles at 
the ultraviolet cutoff scale A. 

Taking into account the non S-wave contributions obtained in appendix A, we obtain the total 
vacuum-energy densities (Tj)in (f7T|l and (T/)out (f7^ - 

<^'>" - ^1 (^) 
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up to the leading order O(A^). 

We rewrite vacuum-energy densities Eq.(IHH) and Eq.(IHSl) as 



,^^v A 5 / Kdr , 

i^tUut ^^^2g-i/2dr6 [ 2TT )■ ^^^^ 

Eg . (18611 indicates that in the shell-volume 47rr^dr, the number of quantum-field states is (Adr/lir) 
and these states carry the energy A. Ea. ()87() shows that the number of quantum-field states 
(Adr/27r) is invariant, the energy of these states receives the gravitational red-shift g^^'^{r) and the 
shell- volume 'iirr'^dr receives the gravitational factor g~^^^{r). 

VII. VACUUM-ENERGY DIFFERENCE AND ITS VARIATION 

We have obtained the vacuum-energy densities (T/)in (|84|) and (T/)out l|85|l . and in this section 
we compute vacuum energies in the presence and absence of gravitational field. With respect to 
a static observer in spacetime 0, whose killing vector and four velocity are and u^, we define 
the vacuum energy as 

£{x) = u^{T^,)u''dll = n''{T^,)dTy, (88) 

and 

^totai = / u^'{T^u)u''d^ = I n^'{T^,)dJ:\ (89) 

where is the hypersurface described by the equation t = Constance, dT,'^ its hypersurface element 
vector, 

dJ:" = n^d^, = f , = 4^, (90) 

and the hypersurface element = ^/hd^x, hij is the spatial metric. 

At the initial time (tin = —St/2), with respect to the static observer O H34|) located in the 
internal flat spacetime 5, the vacuum energy is given by 

fin(r) = 47rr2dr(r/)in, 



6V 2tt 

where the vacuum-energy density {Tl)i^ is approximately given by Ea. H84() . While, at the final time 
(^out = +St/2), with respect to the same static observer O located in the external spacetime 
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S, the vacuum energy is given by 

^out(r) = \^d^x{T})out, 

= A^r\-g„fl^dr{Tl),^,, (92) 

where the vacuum-energy density (T/)out is approximately given by Eq.ljHHJ)- Note that Ea. H92() is 
the same as Eq.(6.184) in the book by Birrell and Davies 

Using Ea. (|M|) for £\n{T) and Ea. (|U5|) for iSout(?')i we approximately obtain the vacuum energy 
£o\\t{i') in the presence of gravitational field, in terms of the vacuum energy fin(r) in the absence 
of gravitational field, 

£ont{r) = g"\r)&,^{r), (94) 

which are the same as the results obtained by using the Heisenberg uncertainty relationship in 
ref. J]. We find that the vacuum energy <?out(?') in the presence of gravitational field is gravitation- 
ally red-shifted from the vacuum energy £ia{r) in the absence of gravitational field. 

Corresponding to the difference (T/(x))difr (j69|l . the difference between the vacuum energy (|93|l 
and vacuum energy (PT|l is 

5£{r) = fout(r) - fi„(r) = {g^'\r) - l)£,^{r) < 0, (95) 

which shows the vacuum energy gets smaller, implying that the vacuum state gains gravitational 
energy when gravitational field is turned on in the time interval 5t = tout — ^in- For r S> 2M, we 
approximately obtain the difference ()95() of vacuum energies, 

5£{r) = fout(r) - fin(r) ^ -^^Ml, (96) 

This indicates an interacting energy due to the vacuum energy £in{r) coupling to the negative 
gravitational potential —M/r. 

For the static observer O located at r in the flat spacetime S, absolute value of the vacuum 
energy £'in(?') (|91|) is not measurable. Analogously, the static observer O located at r in the curved 
spacetime S is not able to measure the absolute value of the vacuum energy <?out('') (ESI)- The 
question is then how to show the difference H95|) of vacuum energies in the presence and absence 
of a static gravitational field on the Earth. In the short letter Q], author suggested to measure the 
Casimir effect at different altitudes above the Earth so as to reveal the difference of vacuum energies 
due to gravitational field acting on the vacuum. Although the gravitational field is static, it varies 
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in the radial position, the static observer can possibly detect the variation of such difference (|96() . 
by measuring the Casimir effect at different altitude ri and r2, 

Casimir Casimir; 

(97) 

where 5iS(r)casimir(?'2) is the difference of vacuum energies in the presence and absence of two 
conducting plates. Such a gravitational effect (jHTJ on the vacuum energy seems too small to be 
seen. 

However, in a gravitational collapsing process approaching to the horizon of a black hole, with 
respect to a static observer located at R{t), where is the collapsing shell and gravitational field 
strongly and rapidly varies, the difference 5£{r) ()95() could be very large and the vacuum state 
gains a large amount of gravitational energy. 



VIII. GRAVITATIONAL COLLAPSING 



In order to be able to analytically study the vacuum-energy gain and in particular vacuum 
decay leading to photon productions in the following sections, we simplify the dynamical process 
of a gravitational collapse and adopt a model with the following approximations: 

• gravitational collapsing of massive and spherical shell that is infinitesimally thin; 

• exactly spherical symmetry in gravitational collapsing; 

• stationary Schwarzschild geometry in the external spacetime S outside the massive shell; 

• flat geometry in the internal spacetime S inside the massive shell. 

This spherical shell has a total mass-energy M and rest mass Mq. In this model, the solution to 
the Einstein equations was studied in refs.j2^ and such a gravitational collapse process can be 
described by the equation, 

6R g{RWK\R)-g{R) 



6t h{R) 

, , ^ 2M 1 

h(R) = r , 

^ ' ART' 



(98) 



where V{t) = 6R/dt, in the unit of c, is the collapsing velocity in the opposite radial direction, and 
the collapsing parameter T = M/Mq. The solution of this equation R = R{t) is the radial location 
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of the collapsing shell at the moment t. At the moment to when the collapsing velocity V{to) = 0, 
the shell, that is at rest and starts to collapse, is located at the radial position: 

, , 1 2M 

Roito) = (99) 

We chose M = WMq and T = 0.0257, Eq.® gives i?o(to) = 10(2M). In Fig.©, we plot the 
collapsing velocity V{t) as a function of R = R/2M, shell's radius in unit of 2M. We find that the 
collapsing process R{t) for Ro{to) > R{t) > 2M undergoes almost in the speed of light, however, 
it becomes very slow when it approaches the horizon: R{t) 2M. The collapsing process takes 
about 0.004 seconds. 

When such collapsing shell R{t) sweeps inwards 6R in the time interval 5t, with respect to 
the static observer located at R{t), the vacuum state changes from |0,in) @ to |0,out) ((TT|) . 
The amplitude of vacuum state transition is given by Ea. H24|) . correspondingly the variations of 
energy-momentum densities are given by Eqs. (|27l(i9)) . The vacuum-energy variation is given by 
Eq.jnSl), 

6£iR) = {g^'\R) - l)£UR) 

= {g^/\R)-l){^^J^5R, (100) 

indicating that vacuum state |0,out) (jllj) gains gravitational energy with respect to |0,in) ©. By 
using the collapsing equation (|98jl , we obtain the rate of vacuum-energy variation per unit of time: 

where and henceforth we take the absolute value of the collapsing velocity V = 6R/6t. This is the 
rate of the vacuum states gaining gravitational energy in gravitational collapsing process. 

In order to see the numbers of this rate in Ea. HlOlf) . we take the ultraviolet cutoff A = Ap and 
convert the natural unit (Ap = 1) into: A^ = ^ = 1.95- lO^^ergs/5.4- 10~^'^sec = 3.6- lO^^egrs/sec, 

^^(R) _ ,f 1 n58 n ,1/2 ^u.-^^R f eg^^s 



4.78-10^«(l-5^/^(i?))— ^ . (102) 
ot ot V sec / 

It is worthwhile to point out that the number 3.6 • lO^^egrs/sec in the rate of the vacuum-energy 
gain is completely determined by natural constants, independently of any free parameters. 

Eas. H102() completely determines the rate of vacuum-energy gain 5£{R)/5t in the spherical shell 
At[R?5R that the collapsing shell sweeps inwards in the time interval St. We plot the rate of 
vacuum-energy variation (gain) 5£{R)/5t in Fig.©. The result shows that the rate 5£{R)/5t 
rapidly increases to 10^"^ erg/sec, as the radius R{t) of the collapsing shell moves, almost in the 
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speed of light, inwards to the horizon. Whereas, in the vicinity of the horizon, the cohapsing 
process becomes slow and the rate 6£{R)/dt decreases and goes to zero. 

The total amount of vacuum-energy gain from gravitational field at the end of gravitational 
collapse is given by integrating Ea. ()lUO() : 

^totai = Jj^(^)=[l)^j^ {9^'\R)-l)5R 

where Ro = 10(2M). The maximum variation of gravitational energy is M/2 in the collapse 
process, which can be derived from differentiating the gravitational potential {—M/r) from r ~ oo 
to r = 2M. 

Due to this vacuum-energy gain 5£{R) (jlOUj) . vacuum states become energetically unstable, 
have to spontaneously undergo a quantum transition to lower energy states via quantum-field 
fluctuations, which leads to particle productions. As a consequence, the vacuum-energy 5£{R) 
H1UU|) gained from gravitational field must be released and deposited in the region from r = 2M to 
r = i?o. 

IX. ENERGY RELEASE AND PHOTON PRODUCTIONS 

Which process of quantum transition releases this vacuum-energy gain 58 {R) H1UU|) ? One of 
possibilities is spontaneous photon emission, analogously to the spontaneous photon emission taking 
place in the atomic physics. In ref. 4], we mentioned the possibility of such a spontaneous photon 
emission can be induced by the four-photon interacting vertex in the Quantum Electromagnetic 
Dynamics (QED). Yet, we have not been able to calculate the rate of such spontaneous photon- 
emission, since we are studying quantum scalar fields, instead of the QED, in curved spacetime. 
In this section, we try to compute the quantum transition amplitude between the initial and final 
vacuum states, corresponding to before and after gravitational field is turned on, at each step of 
gravitational collapsing process. This quantum transition amplitude is related to the probability 
of spontaneous particle ("photon") productions. 

First, we define the invariant scalar product of initial vacuum state and final vacuum 

state </>out(2;) (see Eq.(3.28) in the book by Birrell and Daviesj^): 

(<^out,</'in) = -i I (^ont{x)d^^Ux)d^'' (104) 

<Pontix)df,(j)*^{x) = (l)ontix)df^(j)*^{x) - df^[(j)out{x)](l)*r,{x) 
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where dS'^ = nf^dT,, with a future-directed unit orthogonal to the spacehke hypersurface S and dT, 
is the volume element in S. Since the value of {(j)out, 4'in) is independent of S, we rewrite Eq. (|l()4jl 
as, 

(</'out,(/'in) = -W dT,^4>out{x)dt(t)*n{x) (105) 
J St 

where St is the spacelike hypersurface for t =constance and its element dS* = n*dS is given in 

Eq.jnni). 

With respect to the rest observer O located at the radial position r = R{t) and at the moment 
t, where and when the collapsing shell R{t) sweeps inwards, the gravitational field is turned on for 
the time interval 6t. This time scale 5t is determined by the gravitational collapsing process 5R 
(|98|) . In this time interval 6t, the vacuum state changes from the initial vacuum state (f)in{x)- 

to the final vacuum state (j)ontix)' 

1 „-iu!ntn n . f^W. fn\ ^ ^ D SB. 



e-*-°*°i2.oioWn,mo(^^), r<R-f 



We assume that the time scale of quantum-field transition from the initial vacuum state (j)m{x) 
to the final </>out(3;) induced by the variation of gravitational field is very much shorter than the 
gravitational collapsing time-scale 6t, so that the gravitational field is adiabatically turned on at 
r = R{t), the initial vacuum state (/>in(a^) and final vacuum state <^out(2^) are considered as their 
asymptotic eigenstates (a;o,^0;"io) and {uj,l,m) respectively. 

With these initial state (j)m{x) and final states (j)out{x), Eq. (|l()5j) gives the vacuum to vacuum 
transition amplitudes, when the gravitational field is turned on at r = R{t). These vacuum to 
vacuum transition amplitudes are just the Bogolubov coefficients: 

ajj = (</>out, 0in); Pij = -{(t>out,(t>in), (108) 

where |/3ijP describes the probability of particle productions. Using H106|) and (pox^t <|107|1 . we 
compute the transition amplitude f3ij = {(pout, (ptn) 



W,<Ain) = / dT. {ll>o - UJo)(j)out4>m 

J[0,R-5R/2] 



+ dT, {uj- u;o)(/>out</'in 

J[R-5R/2,R+5R/2] 
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+ / di;*(cj - a;)</)out0m 

J[R+5R/2,oo\ 

dS*(w-a;o)'/'out0in, (109) 

[R-SR/2,R+5R/2] 

where [A, B] indicates the integration zone of in the radial direction, and in the last line of 
equation 



— --'""*"iic.„io(0^o-o(f^)> (110) 



where lo = g^/'^{r)ojQ^ 

Summing over initial states, we obtain the probability of particle productions in final states 
within the energy interval (cj, w + da;), 



R-SR/2 Ag^/^{r) 

where / = Iq and m = mo- To derive Eq. (|llH) we use the orthogonality and closure relations of 
eigenfunctions R^oi^{r)Yig^^{n): 



(1) d^}Yll^{momoi^) = 



mmo ) 



(2) J2 {R.oloir)Ylon.om* R.oloir)yiomom 
= ^6{r - r')6\n - n'); 

(3) / d^'^6{r-r')6\n-n')f{r,n) = f{r',n'). 

Using the rate 6R/5t given by gravitational collapsing equation (jUHJ, we obtain the rate of the 
particle productions. 



R'\R^i{R)\'. (112) 



dtdu; 4.gy^{R) 6t 

This equation gives the rate and spectrum of particle productions, corresponding to the vacuum- 
energy variation 6£{R)/6t (jlUip in the simplest model of gravitational collapsing shell, as described 
in the beginning of section HVIII|) . 
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In order to have an idea of the number of particle creations in a second, we approximately use 
the continuity of functions RujQiiR) and Roji{R), at r = R 

RM-RuoiiR), (113) 
where is infinitesimally thin shell in gravitational collapsing. As a result, we have 



dtduj 4:gy^{R) 6t 

where R^^gi = 2uJoji{uJoR). Using the relation uj = g^^^{R)ujQ and integrating over [0, A] in 
Ea. (|114() for the S-wave (/ = 0), we obtain, 



ot Vsec. / 



where we take A = Ap again. We plot the rate of particle creations dN/dt in terms of ^ = R/2M 
in Fig.©, which shows that the rate 6N/6t rapidly increases to 10^^/sec, as the radius R{t) of the 
collapsing shell moves, almost in the speed of light, inwards to the horizon. Whereas, in the vicinity 
of the horizon, the collapsing process becomes slow and the rate 6N/5t decreases and goes to zero. 
The total number of particles created in the collapse process is given by integrating Ea. (|115|) . 

N = ^{2M)j'\l-gy\R)f6R, 

= OMIApIiMq = 6.69 • 10^^ (116) 

As shown in Figs.© and ©, the rate of vacuum-energy gain, and the rate of particle creations 
are very large, as the collapsing process approaching to the formation of black hole's horizon 
R = 2AI. The total energy output (|1U3() and number H116() of particles created are enormous. 
These qualitatively agree to the characteristics of energetic sources for gamma ray bursts. It is 
indeed interesting that the Planck scale Ap as the ultraviolet cutoff in our proposal neutrally gives 
rise to the characteristics of gamma ray bursts, instead of depending on an arbitrary energy scale. 

The energy of photons spontaneously emitted can be larger than the energy threshold 2me, so 
that electron and positron pairs are produced. These pairs, on the other hand, annihilate into two 
photons. As a consequence, a dense and energetic plasma of photons, electron and position pairs, 

energy and particle- 
number densities of this ^^dyadosphere" can be respectively obtained by Ea. (|lU2jl and Ea. (|115j) . 
see Figs. ()2l3p . The total energy and particle-number of ^^dyadosphere" are given by Ea. H103|) and 
Ea. (|TTHl) . 
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X. COMPARISON WITH SONOLUMINIESCENCE 

Sonoluminiescence is the phenomenon of the intense flashes of Hght emission by the pulsations 
of a gas bubble driven by sound-wave in fluid^^. Such experiments deal with the pulsations of 
bubbles of air in water, driven by a sound wave of frequence of 20-30 KHz. During the expanding 
phase, the bubble radius reaches maximum of order R ~ 4.5^m, followed by a rapid collapse down 
to a minimum radius of order R ~ 0.5/im. The photons are emitted, having a "quasi-thermal" 
spectrum with a "temperature" of several tens of thousands of degrees Kelvin. There are about 
10^ photons emitted per flash, and the time-averaged total power emitted is between 30 and 100 
mW. The photons appear to be emitted a very tiny spatio-temporal region: Estimated flash widths 
vary from less than 35 ps to more than 380ps depending on the gas in the bubble 

The fundamental mechanism of such photon emissions in this phenomenon is still very 
controversial 3] • We do not want to enter these controversial discussions in this article. In this 
section, we attempt to briefly discuss the Schwinger proposal[3| in the connection with our study 
of large photon productions in a gravitational collapse. Schwinger considered this phenomenon of 
photon emissions as the Casimir energy (vacuum energy) E releasing, due to the variation of the 
Casimir energy when a very rapid collapse of dielectric material into a vacuum takes place, 

r drdk 1 , r, / 1 \ , . 

E = - j—r:^-\k\[l-^=\, (117) 



(27r)3 2' ' \^ ' 

where e(r) is dielectric constant and |A:| the vacuum energy of zero-point fluctuation mode. The 
total excess energy is 

1^1 = ^^^•<' (l - ^^j) . (118) 

for a slow varying dielectric constant e, where R is bubble's radius and K \s a. cut-off wavenumber. 
The dielectric constant e ^ 1, with respect to the high-energy modes above the cutoff K. If 
this Casimir energy releasing is completely in form of photon emissions, one identifles the average 
number of photon emissions as, 

[ drdk \ , r- -,\ 

^ = y(^2(^-^)' 

= l-R^K^^e-l). (119) 

The cut-off wavenumber K ~ lO^cm"^ within the ultraviolet region, the energy-budget pi8|) and 
the number of photon emissions N (|119j) agree with the experiments, although the spectrum of 
Sonoluminiescence does not extend to the ultraviolet region. However, Schwinger neither explicitly 
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worked out the mechanism of photon productions nor computed the rate of photon productions in 
the dynamical circumstance that the vacuum-energy variation is very rapid, for very rapid collapse 
of bubble. On the basis of this proposal, there are many further studies in the literatures [13], 
concerning on mechanism of photon productions and other relevant aspects relating the Schwinger 
proposal to the experiment of Sonoluminiescence. We want compare our proposal presented for 
gamma ray bursts in this article with the Schwinger proposal for the phenomenon of Sonolumini- 
escence. Regarding the variation of the vacuum energy, we find that vacuum-energy variation H95|l 
due to turning on gravitational field is similar to vacuum energy variation (jllTj) due to changing 
the dielectric constant. The variation of vacuum energies in both equations is negative (e > 1), 
implying the vacuum state gains energy in both cases. In Ea. ()95() . the vacuum state gains the 
gravitational energy. While, in Ea. (|117|) . the vacuum state gains the sound-wave energy. In both 
cases, the variation of vacuum energy is very rapid, because collapsing processes driven by either 
gravitational field or sound-wave are very rapid. The collapsing velocity ~ c in the gravitational 
collapsing case and R ~ 4March in the Sonoluminiescence case. The cutoff wavenumber is a real 
physical cutoff of its own right that Eqs. (|118lll9|) make physical scenes up to this cutoff. Analo- 
gously, The scale A is a real physical cutoff of its own right that Eas. H94l95|) take into account the 
total variation of vacuum energy, attributed to turning on an external gravitational field. It is not 
an artificial cutoff introduced for regulating calculations of divergent terms, and then removed in 
renormalizable theories. Although Schwinger did not explicitly work out the mechanism of photon 
productions in the phenomenon of Sonoluminiescence, the basic idea for vacuum-energy variation 
and photon productions is similar to that we propose in ref.0| and this article: the vacuum state 
gains (sound- wave/gravitational) energy and becomes unstable and must decay to the lowest energy 
state, releasing the (sound- wave/gravitational) energy it gains. 



XI. CONTRASTION WITH THE HAWKING EFFECT 



It is important to differentiate the Hawking radiation from the effect discussed in this article. 
It is clear that both effects are attributed to an external gravitational field interacting with virtual 
particles in the vacuum. However, they are very different, not only in the phenomena of their 
appearances, but also dynamics of their origins. 

First, we see the aspect of phenomenon. The Hawking radiation is black-body radiation from 
a thermal bath of the temperature l/SvrM determined by the scale M. Particle creations leading 
to the Hawking radiation do not depend on gravitational collapse processes, or in the other words, 
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the Hawking radiation can be created by an external static gravitational field. While, particle 
creations discussed in this article do not have a black-body spectrum and the energy-scale of 
particle creations processes is a ultraviolet cutoff A. Such particle creations strongly depend on the 
gravitational collapsing processes, or in the other words, particle creations discussed in this article 
cannot occur in an external static gravitational field. 

Second, we discuss the aspect of dynamics. There are many elegant scenarios interpreting 
the origin of the Hawking radiation around the horizon of a black hole. In the recent article^], 
based on the context of quantum field theories for particle and antiparticle creations in an external 
gravitational field, author presented a scenario for understanding the origin of quantum radiation 
of the Hawking type when gravitational field is present. Since the general formulation of quantum 
scalar-field theory discussed in sections (JJ) and HHI|) is similar to that in ref.P|, we adopt our 
scenario to contrast two different dynamical origins of the Hawking radiation and particle creations 
discussed in this article. 

As discussed in the ref.Ql, a gravitational field polarizes the vacuum: quantum field fluctua- 
tions (creation and annihilation) of virtual particles and antiparticles (positive and negative energy 
states) in the vacuum are "aligned" by an external gravitational fleld. By this gravitational po- 
larization effect, the vacuum gains gravitational energy. This energy-gain reduces the energy-mass 
gap that is a barrier, preventing virtual particles in the vacuum from tunneling and creating par- 
ticles. As a result, the probability of such quantum tunneling is increase. This causes the vacuum 
decay and creations of particles and antiparticles. This quantum tunneling effect is quantitatively 
described by the imaginary term, i2M/r^ in Eq. (|42)) and imaginary effective action. This is rather 
analogous to the QED vacuum in the presence of an external electric field. The quantum-filed 
fluctuations of charged virtual particles are polarized by the electric field. Particles and antipar- 
ticles are created by the Schwinger mechanism P], when electric fleld strength is strong enough 
to overcome energy-mass gap ~ 2me. The thermal nature of quantum radiation of the Hawking 
type is due to the CTP invariance in the processes of particle creations and annihilations. The 
temperature of thermal radiation is determined by the vacuum-energy gain when gravitational field 
polarizes the vacuum. The reason why the temperature (or vacuum-energy gain) is the order of 
has been discussed in section (|V|) . 

Instead, as discussed in previous sections, the dynamics of particle creations we discussed in 
this article is very similar to the dynamics of mechanism that Schwinger discussed for Sonolumini- 
escence. It is also rather similar to the dynamics of photon-creations of the dynamical Casimir 
effect jisl. For the reasons that gravitational fleld interacting with virtual particles in the vacuum 
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and the variation of gravitational field in the collapse process, the vacuum gains gravitational en- 
ergy and the large variation of vacuum energy occurs in a very short time and small space. Such a 
large vacuum-energy gain makes the vacuum state to be energetically unstable. Via quantum-field 
fluctuations, unstable vacuum state has to transit to lower energetical vacuum-state. Such vacuum 
to vacuum transition releases the gravitational energy, that the vacuum state gains, by spontaneous 
photons emissions. On the basis of its dynamical origin, this process of photon productions clearly 
does not take place in an external static gravitational field, very differently from the thermal radia- 
tion of the Hawking effect. It can be seen from computations in previous sections that this process 
is mainly contributed by the variation of terms —g^^ij^ — g^^k"^. in Eq. (|42|) . during a gravitational 
collapse. This contrasts with the imaginary term i2M/r^ in Ea. (|42l) . describing quantum tunneling 
process for the Hawking effect in a static gravitational field. The energy scales of two processes 
are very different, which have been discussed in section (|V|) . 



XII. SOME REMARKS 



The research of our proposal is at a preliminary step. We adopt the action ^ for scalar fields, 
rather than the vectorial field of electrodynamics dynamics (QED) in curved spacetime. The notion 
of photon productions in the title and text of the present article should be replaced by "photon" 
productions. In appendix A, we make an approximation in computing non S-wave contributions. 
We would like to consider the results (|84I85() as the S-wave contribution only. In addition, we 
adopt an approximate model (infinitesimally thin shell) for describing the process of gravitational 
collapse. Learning a controversy^] on the ultraviolet cutoff introduced in the Schwinger proposal 
for Sonoluminiescence, we need to further strengthen our discussions and arguments on the ultra- 
violet cutoff A and its value in our proposal. As discussed in section > Schwinger introduced the 
ultraviolet cutoff K to agree with the energy budget and particle number of Sonoluminiescence, al- 
though the spectrum of Sonoluminiescence does not extend to this ultraviolet region. Analogously, 
the ultraviolet cutoff A at the Planck scale is introduced in our proposal to be consistent with 
the energy budget and particle number of gamma ray bursts in sections ()VIIHIX|) . although the 
spectrum of gamma ray bursts is O(MeV). From the spectrum of particle creations Eqs. (|112j) or 
H114|) . we find that particle creations are dominated in the low-energy region, since the spectrum is 
approximately related to the function s\v?{ujR) / {loR)"^ for a given value of collapsing radius R. We 
speculate that high-energy particles should lead to multiparticle productions, and total number of 
particles produced is much larger than ()116|) . We are still far from a complete understanding 
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of our proposal for gamma ray bursts. Nevertheless, it is highly deserved to study the proposal 
presented in this article in connection with the origin of gamma ray bursts. 

It is a good analogy to compare our proposal for the origin of gamma ray bursts with Schwinger 
proposal for the origin of Sonoluminiescence. Further experimental and theoretical studies on the 
Schwinger proposal for Sonoluminiescence will definitely help us to have a better understanding of 
the origin of gamma ray bursts in our proposal. Beside, experimental and theoretical studies on 
photon productions in the dynamical Casimir effect are essential for us to further understand the 
mechanism of photon productions in our proposal. 

In literatures [2^. we find that the Schwinger idea for Sonoluminiescence has been applied for 
explaining the origin of gamma ray bursts, on the basis of the variation of dielectric constant during 
a gravitational collapse. These scenarios seem interesting, in particular, in explaining the total 
energy budget of gamma ray bursts. Analogously, conducting electron gas is used as boundary 
conditions for computing the Casimir energy to discuss possible huge output of cosmic energy 
accounting for Quasars [2J|. 

In future work, we expect to be able to study the quantum field theory of electrodynamics 
dynamics (QED) in curved spacetime and use more precise model describing the process of gravita- 
tional collapse, as well as elaborate calculations of vacuum-energy density and rate of gravitational 
energy releasing by spontaneous photon productions. 



I am grateful to R. Ruffini for bring me into this arena of physics and his continuous support. 
I thank to L. Vitagliano for many discussions on the issues of observer and gravitational collapse. 



In this appendix, we calculate the non S-wave (/ / 0) contributions {T^) * in Eq.dlll) and 
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is given by, 





(120) 
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1=1 



K h 



(121) 



where we use the formula (|74p in Ea. (|121() . Using Eas. H74() and (I77() . we compute (T'*)out Eq.ljTSJ, 



^1/2 oo /. 



(2vr) 



/ 9 2M 



r(-f) ^;(; + i) 



(122) 
(123) 



Based on (T/)|^° (tnTj) and (r/)ott for e ^ up to 0(2/e), we find that the terms and 

/(/ + l)/r^ are canceled in the difference {T^^^^ Ea. H69() . 

In the following calculations, we only keep up to the terms that are 0(2/e) in the limit of e ^ 0. 
Thus, considering the difference Ea. l|E^ . we write Ea. (|120() and Ea. H122p as, 



X-'t / out 



1=1 



47rr2 J (27r) ^ 
The summation over "/" in Eqs. (|124j) and (|125j) is given by 

oo 

E(2/ + 1) = 2C(-1)+C(0), 
1=1 

where C('^) is the Riemann zeta function 

oo -j^ 

C(") = E7;r' ^ = 0,1, 2,3,- 



(124) 
(125) 



(126) 



1=1 



(127) 



Resulted from summing over "Z" in Eqs. (|124|) and (|125j) . C(~l) ^^'^ C(0) ^-re divergent. As has been 
discussed, the transverse momenta k± of transverse modes do not directly couple to the gravita- 
tional field, we discard these divergent terms, by assuming these divergent terms are independent 
of gravitational field g{r). To eliminate these divergent terms, we use the reflection formula of 
analytic continued T- and (^-functions for complex variable z = —n + 5 and (5 — > 0, 



n-{)a-z) = ^-^-'''n^)az+i). 



(128) 



As results. We obtain 



/in 



6 ) Anr"^ J J^k) 



33 



1\ 1 /A2^ 



6 / 47rr^ \ 27r / ' 
1\ g /A2\ 



(129) 



47rr2 \^2^; • 

In the second lines of these equations ()129() and ()13U() . we only keep the leading order O(A^), which 
is in accordance with the limit of e ^ 0. 

It should be pointed that the definitions of fc^o (j^n]) and fc^g depend on the angular quantum 
number "/" and this implies that e in Eas. ()121() is /-dependent. The exchanging the order of k^Q- 
integration and /-summation in Eas. H12Ulll21() is not exact for a finite e, and the same problem for 
Eqs. (|122ll23j) On the other hand, the /-dependence in Eq. (|12()|) is dominated by the degeneracy 
term 2/ + 1 and angular momentum term /(/ + 1). Compared with these dominate terms, the 
/-dependence in k^-Q and k'^Q is very weak in terms of the radial wave- function ji^^^ , whose value is 
limited for / oo. Thus, we neglect the /-dependence of (/crOj/Cro) ^^'^ {kr,k'^) in computations. 



We consider this exchanging to be a good approximation for e ^ 0. The S-wave results (|82l83j) are 
free from these approximations. 



[1] S.W. Hawking, Commun. Math. Phys. 43, 199 (1975); G.W. Gibbons and S.W. Hawking, Phys. Rev. D 
15 2752 (1977); and see for examples: 
T. Damour and R. Ruffini, Phys. Rev. D. 15, 332 (1976); 

R. Brout, S. Massar, R. Parentani and Ph. Spindel, Phys. Rep. 260, 329 (1995) references there in; 
G. tHooft, Nucl. Phys. B256 (1985) 727; Int.J.Mod.Phys. All (1996) 4623-4688; 
P. Kraus and F. Wilczek Nucl. Phys. B433 403 (1995), B437 231 (1995); 
M. K. Parikh and F. Wilczek, Phys. Rev. Lett. 85, 5042 (2001) and references therein; 
R. B. Mann, L. Tarasov and A.I. Zelnikov, Class. Quantum Grav. 9 (1992) 1487; 
A. Ghosh and P. Mitra, Phys. Rev. Lett. 73 (1994) 2521; 
D. Kabat and M. J. Strassler, Phys. Lett. B329, (1994) 46; 
L. Susskind and J. Uglum, Phys. Rev. D50 (1994) 2700; 
J.L.F. Barbon and R. Emparan, Phys. Rev. D52 (1995) 4527; 
J.-G. Dcmers, R. Lafrance and R.C. Myers, Phys. Rev. D52 (1995) 2245; 
F. Belgiorno and S. Liberati, Phys. Rev. D53 (1996) 3172; 
V.P. Frolov, D.V. Fursacv and A.L Zelnikov, Phys. Rev. D54 (1996) 2711. 
[2] N.D. Birrell and P.C.W. Davies, "Quantum field theory in curved space" Cambridge University Press 
(1982); 



34 



R.M. Wald "Quantum field theory in curved space and black hole thermodynamics" University of 
Chicago Press (1994); 

V.P. Frolov and I.D. Novikov "Black hole physics" Kluwer Academic Publisher (1998) references there 



[3] S.-S. Xue, |hep-th/0204060[ final lengthy version revised in Oct. 2003. 



[4] S.-S. Xue, Mod. Phys. Lett. A 18, No. 19 (2003) pp. 1325. 

[5] F. Sauter, Zeits. Phys. 69 (1931) 742; W. Heisenberg, H. Euler, Zeits. Phys. 98 (1935) 714; J. Schwinger, 

Phys. Rev. 82 (1951) 664. 
[6] J. Schweppe et al., Phys. Rev. Lett. 51 (1983) 2261; 

I. Ahmad et al., Phys. Rev. Lett. 75 (1995) 2658; 

R. Ganz et al., Phys. Lett. B 389 (1996) 4; 

U. Leinberger et al., Phys. Lett. B 394 (1997) 16; 

R. Bar et al, Nucl. Phys. A 583 (1995) 237; 

S. Heinz et al., Eur. Phys. J. Al (1998) 27; 

D. L. Burke et al., Phys. Rev. Lett. 79 (1997) 1626. 
[7] S. S. Gerstein, Y. B. Zel'dovich, Lett. Nuovo Cimento 1 (1969) 835 and Sov. Phys. JETP 30 (1970) 

358; 

V. S. Popov, T. I. Rozhdestvenskaya, JETP Lett. 14 (1971) 177; 
V. S. Popov, Sov. J. Nucl. Phys. 14 (1972) 257; 
Y. B. Zel'dovich, V. S. Popov, Sov. Phys. Usp. 14 (1972) 673; 
A. Ringwald, Phys. Lett. B 501 (2001) 107; 

R. Alkofer, M. B. Hecht, C. D. Roberts, S. M. Schmidt, and D. V. Vinnik, Phys. Rev. Lett. 87 (2001) 
193902; 

C. D. Roberts, S. M. Schmidt, D. V. Vinnik, Phys. Rev. Lett. 89 (2002) 153901; 

K. Kajantie, T. Matsui, Phys. Lett. B 164 (1985) 373; 

G. Gatoff, A. K. Kerman, T. Matsui, Phys. Rev. D 36 (1987) 114; 



S. A. Smolyansky et al., |hep-ph/9712377| 

Y. Kluger, J. M. Eisenberg, B. Svetitsky, F. Cooper, E. Mottola, Phys. Rev. Lett. 67 (1991) 2427; 
Y. Kluger, J. M. Eisenberg, B. Svetitsky, F. Cooper, E. Mottola, Phys. Rev. D 45 (1992) 4659; 
F. Cooper, J. M. Eisenberg, Y. Kluger, E. Mottola, B. Svetitsky, Phys. Rev. D 48 (1993) 190; 
J. C. R. Bloch et al., Phys. Rev. D 60 (1999) 116011; 



D. V. Vinnik et al., Eur. Phys. J. C 22 (2001) 341, |nucl-th/0103073 
[8] T. Damour, R. Ruflini, Phys. Rev. Lett. 35 (1975) 463; 

G. Preparata, R. Ruffini, S.-S. Xue, A&A 338 (1998) L87 and J.Korean Phys.Soc. 42 (2003) S99-S104; 

R. Ruflini, L. Vitaghano and S.-S. Xue Phys. Lett. B559 (2003) 12-19 and B (2003). 
[9] S.-S. Xue, Phys. Rev. D68. (2003), 013004. 
[10] G.L. Israel, T. Oosterbroek, L. Stella, S. Campana, S. Mereghetti and A.N. Parmar, Astro- 



35 



phys. J. Lett. 560, L65 (2001). 
[11] H.B.G. Casimir, Proc. Kon. Ned. Akad. Wetenschap., ser. B, vol. 51, p. 793, 1948. see also M. Fierz, 

Helv. Phys. Acta, vol. 33, p. 855, 1960. 
[12] M.J. Sparnnaay, Physica, vol. 24, p. 751, 1958; 

S.K. Lamoreaux, Phys. Rev. Lett. 78 (1957) 5; U. Mohideen and A. Roy, ibid 81 (1998) 4549. 
[13] see for examples: 

Astrid Lambrecht, Marc- Thierry Jaekel and Serge Reynaud, Phys. Rev. Lett. 77 (1996), 615; 
Ramin Golestanian and Mehran Kardar, Phys. Rev. A 58 (1998), pp. 1713-1722; 
Diego A. R. Dalvit and Francisco D. MazziteUi, Phys. Rev. A 57 (1998), pp. 2113-2119; 
Diego A. R. Dalvit and P. A. Maia Neto, Phys. Rev. Lett. 84 (2000), 798-801; 
P. A. Maia Neto and D. A. R. Dalvit, Phys. Rev. A 62 (2000), 042103 (11 pp.); 
G. Plunien, R. Schtzhold, and G. Soff, Phys. Rev. Lett. 84 (2000), 798-801882-1885. 
[14] B.P. Barber, R.A. HiUer, R. Lofstedt, S.J. Putterman and K. Weniger Phys. Rep. 281, 65-143 (1997); 
B. Gompf, R. Gunther, G. Nick, R. Pecha, and W. Eisenmenger, Phys. Rev. Lett. 78, 1405 (1997); 
R.A. Hiller, S.J. Putterman, and K.R. Weninger, Phys. Rev. Lett. 80, 1090 (1998); 

D. Lohse and S. Hilgenfeldt, J. Chem. Phys. 107 (1997), 6986; 
T.J. Matula and L.A. Crum, Phys. Rev. Lett. 80, 865 (1998); 

J. A. Ketterling and R.E. Apfel Phys. Rev. Lett. 81, 4991 (1998); 

J.B. Young, T. Schmiedel, Woowon Kang Phys. Rev. Lett. 77, 4816 (1996); 

B. A. DiDonna, T.A. Witten, J.B. Young, Physica A 258, 263 (1998); 

S. Hilgenfeldt, D. Lohse, and W.C. Moss, Phys. Rev. Lett. 80, 1332 (19980). 
[15] see for examples: 

E. SassaroH, Y. N. Srivastava, and A. Widom, Phys. Rev. A 50, 1027 (1994); 

C. Eberlein, Phys. Rev. A 53, 2772 (1996), and Phys. Rev. Lett. 76, 3842 (1996) and 77, 4961 (1996); 
78, 2269 (1997); 

C. E. Carlson, C. Molina-Pan's, J. Perez-Mercader, and M. Visser, Phys. Lett. B395 , 76 (1997); Phys. 
Rev. D 56, 1262 (1997); 

C. Molina-Paris and M. Visser, Phys. Rev. D 56, 6629 (1997); 
C.S. Unnikrishnan and S. Mukhopadhyay, Phys. Rev. Lett. 77, 4960 (1996); 
A. Lambrecht, M-T. Jaekel, and S. Reynaud, Phys. Rev. Lett. 78, 2267 (1997); 
N. Garcia and A.P. Levanyuk, Phys. Rev. Lett. 78, 2268 (1997); 

M. Visser, S. Liberati, F. Belgiorno, and D.W. Sciama, Physical Review Letters 83 (1999) 678-681; 

S. Liberati, M. Visser, F. Belgiorno, and D.W. Sciama, Journal of Physics A33 (2000) 2251-2272; 

S. Liberati, M. Visser, F. Belgiorno, and D.W. Sciama, Physical Review D61 (2000) 085023 ; 

S. Liberati, M. Visser, F. Belgiorno, and D.W. Sciama, Physical Review D61 (2000) 085024; 

C. E. Carlson, C. MoUna-Paris, J. Perez-Mercader, and M. Visser, Phys. Lett. B 395, 76-82 (1997), 

C. E. Carlson, C. Mohna-Paris, J. Perez-Mercader, and M. Visser, Phys. Rev. D56, 1262 (1997), 



36 



C. Molina-Paris and M. Visser, Phys. Rev. D56, 6629 (1997), 

J. Holzfuss, M. Ruggeberg, and A. Billo, Phys. Rev. Lett. 81, 5434 (1998), 

K. Milton, Quantum field theory under the influence of external conditions, edited by M. Bordag, 

(Tuebner Verlagsgesellschaft, Stuttgart, 1996), pages 13-23; 

K. Milton and J. Ng, Phys. Rev. E55 (1997) 4207 and E57 (1998) 5504; 

[16] I thank to I. Swason, who brought the Schwinger proposal for Sonoluminiescence|22j into my knowledge 

after he read my short letter jJl that briefly presented the idea of this proposal. This enables me to learn 

the phenomenon of Sonoluminiescence and relevant literatures on the issue. This encourages me to 

write the present lengthy article for detailed discussions on our idea. 
[17] M. H.P.M. van Puttcn, Phys. Rev. Lett. 84, 3752 (2000); M. H.P.M. van Putten and Amir Levinson, 

Science 295 1874 (2002). 
[18] S. Koide, K. Shibata, T. Kudoh and D. L. Meier, Science 295 1688 (2002). 
[19] T. Piran, 1999, Phys. Rep. 314, 575; 

P. Meszaros, Annu. Rev. Astron. Astrphys. 2002. 40:137-69. 
[20] R. Ruffini, C.L. Bianco, S.-S. Xue, P. Chardonnet and F. Fraschetti, International Journal of Modern 

Physics D, Vol 12, (2003) 173. 
[21] W. Israel, // Nuovo Cimento 44B (1967) 1-14; 

V. De la Cruz, W. Israel, // Nuovo Cimento 51A (1967) 744-759; 

R. Ruflini, L. Vitagliano, Phys. Lett. B 545, (2002) 233, and Int. J.Mod.Phys. D12 (2003) 121-128; 
Christian Cherubini, Remo Ruflini, Luca Vitaghano, Phys.Lett. B545 (2002) 226-232. 
[22] J. Schwinger, Proc. Nat. Acad. Sci. 89, 4091-4093 (1992); 
J. Schwinger, Proc. Nat. Acad. Sci. 89, 11118-11120 (1992); 
J. Schwinger, Proc. Nat. Acad. Sci. 90, 958-959 (1993); 
J. Schwinger, Proc. Nat. Acad. Sci. 90, 2105-2106 (1993); 
J. Schwinger, Proc. Nat. Acad. Sci. 90, 4505-4507 (1993); 
J. Schwinger, Proc. Nat. Acad. Sci. 90, 7285-7287 (1993); 
J. Schwinger, Proc. Nat. Acad. Sci. 91, 6473-6475 (1994); 



[23] C. Carlson, T. Goldman, and J. Perez-Mercader, EuroPhys. Lett. 36 (1996) 637; 

C. Carlson and I. Swanson, Foundations of Physics, Vol. 30, No. 5, 2000. 
[24] I.Yu Sokolov, Phys. Lett.A, v.223, p.l63, 1996. 



37 




R/2M 



FIG. 1: The velocity of gravitational collapsing shell in the unit of c as a function of radius R in the unit of 
2M. 
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FIG. 2: The rate of vacuum-energy gain (ergs/sec) as a function of radius R in the unit of 2M. 
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FIG. 3: The rate of particle creations (1/sec) as a function of radius R in the unit of 2M. 



